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ABSTRACT 

Results  from  a  wind  tunnel  test  of  two  delta  wings 
in  close  proximity  are  presented  and  compared  with 
predictions  from  a  vortex  lattice  method.  The  wake 
effects  on  lift  are  slightly  overpredicted  when  the 
aircraft  overlap  in  the  spanwise  direction.  The 
wake  induced  effects  on  pitching  and  rolling 
moment  are  well  predicted.  A  maximum  induced 
drag  reduction  of  25%  is  measured  on  the  trail 
aircraft,  compared  with  a  40%  predicted  reduction. 

NOMENCLATURE 

A  Aspect  ratio 

b  Wing  span 

Ca  Axial  force  coefficient 

Cl  Lift  coefficient 

Cd  Drag  coefficient 

CDo  Zero  lift  drag  coefficient 

Coi  Induced  drag  coefficient 

Ci  Rolling  moment  coefficient 

Cm  Pitching  moment  coefficient 

Cn  Normal  force  coefficient 

Kjj  Induced  drag  factor 

x  Longitudinal  distance,  Fig  2 

y  Lateral  distance,  Fig  2 

z  Vertical  distance,  Fig  2 

a  Angle  of  attack 

*  Aerospace  Engineer,  Associate  Fellow  AIAA 

**  Aerospace  Engineer,  Senior  Member  AIAA 
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INTRODUCTION 

Many  military  missions  involve  formation  flight.  In 
some  instances,  there  is  significant  aerodynamic 
coupling  between  aircraft  arising  from  the 
proximity  of  one  aircraft  to  wakes  generated  by 
others  in  the  formation.  While  aerial  refueling  is 
the  most  common  example,  a  new  one  is  currently 
under  study,  formation  flight  for  drag  reduction. 
This  is  now  considered  practical  due  to  the  ability 
to  determine  aircraft  location  with  high  precision 
and  advances  in  automatic  control  theory.  If  the 
aircraft  are  appropriately  positioned,  upwash  from 
wing  tip  vortices  can  reduce  the  drag  of 
neighboring  aircraft,  increasing  range  or 
endurance1.  It  is  generally  accepted  that  this  drag 

reduction  is  one  of  the  reasons  that  many  migratory 

^  2 

bird  species  fly  in  closely  spaced  flocks”. 
Photographic  studies  of  Canadian  Geese  have 
found  that  the  average  lateral  spacing  between 

adjacent  birds  is  very  close  to  the  optimum 

3 

predicted  by  simple  aerodynamic  theory  . 

Most  of  the  research  in  this  area  has  concentrated 

on  analytical  modeling  using  potential  flow 

4  5 

techniques  ’  .  Very  few  experimental  results  are 
reported  in  the  open  literature.  Brown  et  al6 
conducted  tow  tank  tests  using  a  3%  scale  Boeing 
747  model  as  the  wake  generator  and  a  generic 
fighter  as  the  trail  aircraft.  The  trail  to  lead  aircraft 
wing  span  ratio  was  0.20.  Forces  were  measured 
on  the  trail  aircraft  only,  with  lift-to-drag  ratio 
increases  of  up  to  50%  measured  at  some 
conditions.  Beukenberg  and  Hummel  conducted 
two  flight  tests  in  an  attempt  to  demonstrate  a 
power  saving  under  realistic  conditions.  In  their 
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second  test,  a  Dornier-228  acted  as  a  vortex 
generator  for  a  following  Dornier-28.  Using  an 
automatic  control  system  to  maintain  proper 
position,  an  average  power  reduction  of  10%  was 
measured  on  the  trail  aircraft  over  a  150  sec. 
interval. 

There  have  been  many  wind  tunnel  tests  studying 
the  impact  of  wake  vortices  on  aircraft  far 
downstream.  Most  of  these  center  on  the  study  of 
wake  vortex  hazards  for  commercial  aircraft  and  are 
conducted  in  large  facilities  such  as  the  80x120 
tunnel  at  the  NASA  Ames  Research  Center  or  the 
30x60  facility  at  the  Langley  Research  Center9. 
Bloy  et  al10  investigated  formations  typical  of  aerial 
refueling  with  the  receiver  less  than  one  wing  span 
downstream  of  the  tanker.  Their  tests  were 
conducted  in  a  small  wind  tunnel  with  a  relatively 
large  wing  span  to  tunnel  width  ratio  of  0.7. 
Significant  wall  interference  effects  were  uncovered 
by  taking  measurements  in  both  open  and  closed 
test  sections. 

The  flow  survey  rig  within  the  Langley  30x60  Full 
Scale  Wind  Tunnel  has  recently  been  modified  to 
accommodate  an  aircraft  model  with  an  internal 
force  balance  al.  The  present  work  discusses  test 
results  from  a  formation  of  two  closely  spaced 
tailless  aircraft  using  this  rig  and  compares  these 
with  predictions  from  a  vortex  lattice  method. 

DISCUSSION 
WIND  TUNNEL  SET-UP 

The  wind  tunnel  tests  were  conducted  in  the 
Langley  Full  Scale  Wind  Tunnel,  currently  operated 
by  Old  Dominion  University.  The  test 
configuration  was  a  Lockheed  tailless  aircraft 
consisting  of  a  65  degree  delta  wing  with  a 
sawtooth  trailing  edge  with  sweep  angles  of  y  25 
degrees.  It  is  a  single  engine  design  with  a  two 
narrow  inlets  on  the  lower  surface.  The  inlets  were 
blocked  for  the  present  test. 

Two  1/13  scale  models  were  tested.  These  were 
manufactured  several  years  apart  by  different  firms, 
so  slight  differences  are  present.  They  will  be 
referred  to  as  models  “A”  and  “B”.  Model  B  was 
manufactured  later  and  has  several  bumps  and 
exposed  rods  to  accommodate  internally  actuated 


control  surfaces.  Model  B  was  used  as  the  trail 
model  during  all  of  the  formation  tests,  and  was 
mounted  using  an  aft  sting  which  was  connected  to 
a  vertical  rod.  This  rod  was  connected  to  the  wind 
tunnel’s  primary  model  support  system.  Model  A 
was  mounted  from  the  top  on  a  newly  built 
apparatus  which  was  attached  to  the  existing  flow 
survey  carriage.  The  apparatus  consists  of  an  arm 
which  projects  forward  of  the  survey  carriage, 
under  which  a  vertical  sting  and  angle  of  attack 
control  rod  are  attached.  This  arrangement  was 
selected  in  an  attempt  to  minimize  rig  interference 
effects  on  model  B.  A  photograph  of  the 
experimental  apparatus  is  shown  in  Figure  1 . 


FIGURE  1.  Wind  Tunnel  Models  on  Apparatus 


The  relative  height  of  the  two  models  was  achieved 
by  varying  the  position  of  Model  B.  This  was 
accomplished  by  moving  the  sting  along  the  vertical 
rod  which  was  attached  to  the  primary  support 
system.  Relative  lateral  and  longitudinal  position 
were  achieved  by  moving  the  flow  survey  rig  to 
which  Model  A  was  attached.  All  runs  were 
conducted  at  a  tunnel  dynamic  pressure  of  5  psf, 
which  corresponds  to  a  speed  of  approximately  65 
ft/sec. 

Both  models  were  equipped  with  internal  six 
component  strain  gauge  balances.  The  estimated 
accuracy  of  the  balance  measurements,  based  on 
pretest  balance  calibrations  expressed  in  non- 
dimensional  form  using  the  reference  area  and 
length  of  the  study  vehicle  is  presented  in  Table  1. 
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Additional  details  on  the  experimental  set-up  and 
data  reduction  are  given  by  Gingras1  . 

PREDICTION  METHOD 

12 

The  planar  vortex  lattice  code  HASC95  was  used 
to  predict  the  wake  induced  effects  on  the  trail 
aircraft.  Each  aircraft  was  modeled  with  540 
panels,  using  36  in  the  spanwise  direction  and  15  in 
the  chordwise  direction.  The  spanwise  panels  were 
evenly  spaced,  cosine  spacing  was  used  for  the 
chordwise  elements.  Care  was  taken  to  ensure 
correct  alignment  of  vortex  filaments  and  control 
points  when  the  aircraft  overlapped  in  the  lateral 
direction.  HASC95  uses  the  flat  wake 
approximation,  i.e.,  the  trailing  vortex  legs  extend 
to  infinity  parallel  to  the  x-y  plane,  and  do  not  vary 
with  angle  of  attack.  One  result  of  this  is  that  in  the 
linear  lift  region,  the  wake  induced  lift,  rolling 
moment  and  pitching  moment  increments  are 
invariant  with  the  angle  of  attack  of  the  trail 
aircraft.  The  HASC95  representation  of  the 
vehicles  is  shown  in  Figure  2. 


Lead  aircraft 


FIGURE  2.  Vortex  Lattice  Representation  of  Vehicles 
13 

Rossow  has  shown  that  the  vortex  lattice  method 
is  extremely  accurate  for  predicting  wake  induced 
effects  when  the  span  of  the  downstream  wing  is 
substantially  less  than  that  of  the  wake  generating 
wing.  When  the  spans  of  the  two  wings  are 
comparable,  as  in  the  present  tests,  accuracy  was 
degraded  somewhat  but  trends  were  still  very  well 
predicted.  The  reduction  in  accuracy  was  attributed 

to  distortions  of  the  wake  cause  by  the  presence  of 

13 

the  downstream  wing.  Bloy  et  al  also  found  good 
agreement  with  vortex  lattice  predictions  when 
compared  to  their  open  test  section  results. 

ISOLATED  MODEL  RESULTS 


The  lift  coefficient  of  models  A  and  B  in  isolation 
is  shown  in  Figure  3,  along  with  the  HASC95 
prediction.  For  the  formation  flight  tests,  model  A 
used  a  top  mounted  sting  and  model  B  used  an  aft 
mount.  Model  A  was  tested  in  isolation  using  both 
aft  and  top  mounts  to  assess  sting  effects  and  allow 
comparison  with  prior  tests  which  used  an  aft 
mounted  sting.  Also  shown  are  data  from  two  prior 
tests  1415  of  model  A  using  an  aft  mounted  sting. 
This  was  the  first  test  with  model  B.  The 
differences  between  the  models  are  at  least  as  large 
as  any  facility  or  mounting  effects.  Model  B  on  the 
aft  mount  has  slightly  higher  lift  than  any  of  the 
model  A  results.  The  HASC95  prediction  is  very 
good  below  10  degrees  angle  of  attack  and  is 
slightly  low  at  higher  angles  of  attack. 


Drag  characteristics  of  the  models  in  isolation, 
along  with  HASC95  predictions,  are  shown  in 
Figure  4.  The  two  models  show  very  similar 
results,  with  Model  A  having  a  slightly  larger  zero 
lift  drag  than  model  B.  The  experimentally 
measured  zero  lift  drag  of  model  B  was  added  to 
the  HASC95  prediction.  Three  predictions  are 
shown,  corresponding  to  three  levels  of  assumed 
leading  edge  suction,  or  flow  attachment.  Full 
leading  edge  suction  corresponds  to  a  fully  attached 
flow  condition,  with: 

Cqi  =CL2  / nA 

Zero  leading  edge  suction  corresponds  to  a  fully 
separated  flow  condition,  with: 
cDi  =  CL  tan« 

At  low  lift  coefficients,  the  drag  compares  well  with 
the  full  suction  prediction.  As  lift  increases,  some 
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flow  separation  is  evident  due  to  the  rise  in  drag 
above  the  full  suction  prediction.  The  drag  never 
reaches  the  value  corresponding  to  fully  separated 
flow  case.  Several  values  of  suction  were  run  in  an 
attempt  to  match  the  value  measured  at  8  degrees 
angle  of  attack,  which  is  the  reference  case  for  most 
of  the  formation  flight  effects.  A  value  of  60% 
attained  suction  gives  a  very  good  match,  as  shown. 


FIGURE  4.  Drag  Polars  of  Models  in  Isolation 


FORMATION  FLIGHT  RESULTS 

LATERAL  SPACING  -  LIFT  AND  DRAG 

The  effect  of  lateral  spacing  on  the  wake  induced 
lift  is  shown  in  Figure  5.  These  results  are  for  a 
longitudinal  spacing  of  two  spans  and  zero  vertical 
spacing  (aircraft  in  the  same  plane  relative  to  the 
oncoming  wind).  The  lead  aircraft  is  at  eight 
degrees  angle  of  attack  in  all  cases.  Results  are 
shown  for  the  trail  aircraft  at  zero  and  eight  degrees 
angle  of  attack.  Maximum  lift  loss  occurs  with  the 
aircraft  directly  in  trail  when  it  is  completely 
enveloped  by  the  downwash  from  the  lead.  As  the 
trail  aircraft  moves  outboard,  the  lift  loss  reduces 
and  upwash  is  encountered  which  reduces  once  the 
trail  aircraft  has  passed  the  lead’s  wing  tip  (y/b>l). 

The  HASC95  prediction  shows  the  correct  trend, 
and  the  magnitude  of  the  lift  is  well  predicted  for 
lateral  spacings  greater  than  2/3  span.  The  lift  loss 
is  underpredicted  for  smaller  spacings.  This  may  be 
due  to  a  reduction  in  local  dynamic  pressure 
immediately  behind  the  lead  aircraft,  which 
HASC95  ignores.  The  wake  induced  lift  is  not 


affected  by  the  angle  of  attack  of  the  trail  vehicle. 
This  is  consistent  with  the  HASC95  results. 

Lift-to-drag  ratio  as  a  function  of  lift  for  several 
lateral  spacings  is  shown  in  Figure  6.  The  aircraft 
in  isolation  is  shown  as  the  solid  line,  with 
formation  flight  results  shown  as  dashed  lines. 
These  data  were  taken  from  angle  of  attack  sweeps 
of  the  trail  aircraft  with  the  lead  at  a  fixed  angle  of 
attack.  Drag  reduction  comparisons  are  most 
meaningful  for  the  case  where  both  aircraft  are  at 
the  same  lift  coefficient.  Here,  a  drag  reduction  of 
almost  15%  is  attained  for  the  trail  aircraft. 


FIGURE  5.  Effect  of  lateral  spacing  on  wake  induced 
lift,  x/b=2,  z/b=0. 


FIGURE  6.  Formation  flight  lift-to-drag 
characteristics,  x/b=2,  z/b=0. 

In  the  HASC95  runs,  both  aircraft  are  at  the  same 
angle  of  attack,  so  the  lift  of  the  two  vehicles  is 
almost  always  different.  Prandtl’s  method16  was 
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used  to  extract  the  induced  drag  effect  for  cases  of 
differing  lift  on  the  lead  and  trail  wings: 

CD2  =  cD,o  +  KocL2 2  +  Kl2cL,lcL2 
Here,  K„  represents  the  self  induced  drag  and  K12 
represent  the  effect  of  wing  1  on  wing  2.  The 
induced  drag  change  due  to  formation  effects  can 
be  easily  calculated: 

K\2  _  CL2  cD2~cD,o  | 

ko  cLAy  k0cL22 

The  same  procedure  was  used  to  extract  the 
induced  drag  reduction  from  the  experimental  data. 
Kq  was  determined  from  a  curve  fit  of  the  isolated 
aircraft  lift  and  drag  results.  If  two  untapered  wings 
are  physically  joined,  the  induced  drag  change  of 
the  added  wing  (K12/K0)  is  -1  or  100%  reduction. 
The  HASC95  results  for  the  fully  attached  flow 
case  (dashed  curve)  indicate  a  maximum  reduction 
of  about  70%  at  a  spacing  of  90%  span,  indicating  a 
slight  overlap  of  the  wings.  Simple  horseshoe 
vortex  theory  predicts  a  maximum  reduction  of 
100%  at  7t/4  (78%)  span.  The  reduction  is  smaller 
for  the  present  case  due  to  the  high  degree  of  wing 
taper.  The  HASC95  results  for  60%  attained 
suction  show  a  predicted  drag  reduction  of  40%, 
also  occurring  at  90%  span  spacing.  The  test  data 
show  a  maximum  reduction  of  about  25%, 
occurring  between  80  and  90  %  span  spacing. 
Although  the  overall  trend  of  the  data  matches  the 
prediction,  the  reduction  is  about  15%  less  than 
predicted  across  the  board.  One  possible 
explanation  for  this  discrepancy  is  that  the  induced 
velocities  at  the  trail  aircraft  wing  tip  are  high 
enough  to  significantly  increase  flow  separation. 
Using  a  pair  of  fully  rolled  up  vortices  to  represent 
the  wake,  the  estimated  induced  upwash  angle  at 
the  tip  of  the  trail  aircraft  is  8  degrees  for  the 
present  test  conditions.  Added  to  the  geometric 
angle  of  attack  of  4-10  degrees,  this  almost 
certainly  increases  separation. 


Lateral  Spacing,  Ay/b  (leader  to  left) 
FIGURE  7.  Effect  of  lateral  spacing  on  induced  drag 
reduction,  x/b=2,  z/b=0. 

LATERAL  SPACING  -  PITCH  AND  ROLL 

The  effect  of  lateral  spacing  on  the  wake  induced 
pitching  moment  is  shown  in  Figure  8.  For 
reference,  a  maximum  (30  degree)  elevon 
deflection  gives  a  pitching  moment  increment  of  - 
0.040,  so  the  maximum  induced  moment 
corresponds  to  about  a  50%  control  deflection. 
There  are  three  distinct  peaks  in  the  incremental 
moment,  occurring  at  approximately  one-quarter 
span,  two-thirds  span,  and  wingtip  spacing. 
HASC95  does  an  excellent  job  of  predicting  the 
size  and  location  of  these  peaks. 

The  pitching  moment  data  in  figure  8  show  a  small 
effect  of  angle  of  attack  when  the  wings  overlap 
(y/b<  1 ),  with  an  increased  nose  up  moment  with 
angle  of  attack.  The  zero  angle  of  attack  data  is 
more  accurately  predicted  by  HASC95.  If  it  is 
assumed  that  the  wake  from  the  lead  is  parallel  to 
the  velocity  vector,  then  at  zero  angle  of  attack  the 
trail  aircraft  is  perfectly  aligned  with  the  lead’s 
wake,  which  is  the  HASC95  representation. 
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FIGURE  8.  Effect  of  lateral  spacing  on  wake  induced 
pitching  moment,  x/b=2,  z/b=0. 

To  examine  the  changes  in  sign  of  the  pitching 
moment  result,  pressure  distributions  were  obtained 
from  HASC95.  These  are  shown  for  lateral 
spacings  of  1/3  and  2/3  span  in  Figure  9  and  10. 
They  were  obtained  by  subtracting  the  pressures  of 
an  isolated  vehicle  at  eight  degrees  angle  of  attack 
from  the  pressures  of  both  vehicles  while  in 
formation.  Dashed  curves  indicate  regions  of 
downwash  while  solid  lines  represent  upwash.  The 
contours  are  shown  at  increments  of  0.02  in  terms 
of  Cp.  Although  there  is  a  small  effect  of  the  trail 
aircraft  on  the  lead,  it  cannot  be  seen  at  this 
resolution  in  Cp.  The  black  dots  indicate  the 
moment  reference  center. 

Figure  9  shows  a  downwash  region  that  is  nearly 
aligned  with  the  lead’s  right  wing  tip.  The  upwash 
is  mostly  aft  of  the  moment  center,  with  a  large 
downwash  at  the  nose.  This  results  in  the  nose 
down  moment  shown  in  Figure  8.  Figure  10  shows 
a  small  downwash  region  at  the  left  wing  tip  of  the 
trail,  with  only  upwash  forward  of  the  moment 
reference  center.  This  results  in  the  maximum  nose 
up  moment  peak  shown  in  Figure  8.  It  is  interesting 
to  note  that  the  approximate  dividing  line  between 
downwash  and  upwash  for  the  2/3  span  case  is  not 
aligned  with  the  lead’s  wingtip,  but  with  a  point 
further  inboard. 


FIGURE  9.  Wake  induced  pressure  contours  at 
y/b=  1/3,  x/b=2,  z/b=0. 


FIGURE  10.  Wake  induced  pressure  contours  at 
y/b=2/3,  x/b=2,  z/b=0. 


The  effect  of  lateral  spacing  on  the  wake  induced 
rolling  moment  lift  is  shown  in  Figure  11.  For 
reference,  a  maximum  (30  degree)  elevon 
deflection  gives  a  pitching  moment  increment  of  - 
0.022,  so  the  maximum  induced  moment 
corresponds  to  about  a  100%  control  deflection. 


FIGURE  1 1 .  Effect  of  lateral  spacing  on  wake 
induced  rolling  moment,  x/b=2,  z/b=0. 

VERTICAL  SPACING 


The  effect  of  vertical  spacing  on  the  wake  induced 
lift  is  shown  in  Figure  12.  Results  are  shown  for  a 
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longitudinal  spacing  of  two  spans  and  three  lateral 
spacings,  0.25,  0.75  and  1.25  spans.  Both  the  lead 
and  trail  aircraft  are  at  eight  degrees  angle  of  attack. 
The  scale  of  the  vertical  axis  matches  that  shown 
for  the  lateral  spacing  comparisons.  The  trends  are 
very  well  predicted,  while  the  magnitudes  are  not. 
In  all  cases,  the  maximum  lift  changes  correspond 
to  the  in-plane  (Az/b=0).  For  some  of  the  0.75  span 
results,  the  same  position  was  tested  during 
multiple  runs.  These  additional  points  are  very 
close  to  one  another,  indicating  good  repeatability 
of  the  induced  effects. 


FIGURE  12.  Effect  of  vertical  spacing  on  wake 
induced  lift,  x/b=2. 

The  effect  of  vertical  spacing  on  the  wake  induced 
pitching  and  rolling  moments  are  shown  in  Figure 
13-14.  For  the  most  part,  as  with  the  lift,  the 
maximum  changes  occur  for  the  in-plane  case.  The 
pitching  moment  results  are  very  well  predicted  in 
all  cases.  The  rolling  moment  results  are  well 
predicted  for  the  0.25  and  1.25  span  spacing  cases. 
For  the  0.75  span  case,  the  results  are  mixed.  The 
data  show  a  continually  increasing  moment  with 
height,  while  the  prediction  shows  virtually  no 
effect  of  height. 


FIGURE  13.  Effect  of  vertical  spacing  on  wake 
induced  pitching  moment,  x/b=2. 


FIGURE  14.  Effect  of  vertical  spacing  on  wake 
induced  pitching  moment,  x/b=2. 


POSITIONAL  STABILITY 

One  important  aspect  of  formation  flight  is 
positional  stability.  Positive  stability  results  if  the 
wake  induced  forces  tend  to  return  the  trail  aircraft 
to  its  original  position  when  disturbed.  About  a 
dozen  positional  stability  derivatives  exist,  many 
have  been  studied  by  Bloy  et  al10.  For  the  current 
tests,  three  will  be  examined,  change  in  lift  and 
pitching  moment  with  height,  and  change  in  roll 
with  lateral  separation.  For  the  sign  convention 
selected,  these  are  stable  for: 
dCi  /  dz  <  0 
dCm  /dz  <  0 
dCj  /  dy  <  0 

Derivatives  involving  yawing  moment  and  side 
force  were  effectively  zero  due  to  the  tailless 


7 

AMERICAN  INSTITUTE  OF  AERONAUTICS  AND  ASTRONAUTICS 


configuration.  Apparatus  limitations  prevented 
study  of  other  derivatives  such  as  change  in  roll 
with  bank  angle. 

Stability  boundaries  for  the  change  in  lift  with 
height  derived  from  HASC95  results  are  shown  in 
Figure  15.  Boundaries  extracted  from  the 
experimental  data  are  shown  in  Figure  16.  The 
FIASC95  results  are  symmetric  with  respect  to  both 
the  horizontal  and  vertical  planes  of  symmetry  of 
the  lead  aircraft.  Sufficient  data  were  only  taken  on 
one  side  of  the  aircraft  to  generate  boundaries.  As  a 
result,  results  are  only  shown  for  the  trail  aircraft  to 
the  right  of  the  lead.  The  experimentally  measured 
boundaries  should  be  symmetric  with  respect  to  the 
vertical  plane  of  symmetry. 


a  small  stable  region  beneath  the  vertical  plane  of 
symmetry  outboard  of  wingtip  (1.0)  spacing. 

Stability  boundaries  for  the  change  in  pitching 
moment  with  height  are  shown  in  Figures  17  and 
18.  FIASC95  predicts  six  distinct  stability  regions, 
which  are  separated  by  two  4-corner  boundaries, 
located  at  0.40  and  0.80  span  spacing.  The  data 
show  the  same  general  trends  as  the  prediction.  As 
with  the  lift  prediction,  the  data  do  not  show  4- 
corner  boundaries.  Instead,  the  region  of  instability 
is  continuous  with  channels  separating  the  stable 
regions.  In  addition,  the  data  show  a  small  region 
of  instability  directly  behind  and  below  the  aircraft 
that  was  not  predicted.  The  Appendix  gives  a 
summary  of  the  longitudinal  dynamic  modes  of  the 
aircraft  in  formation  flight. 
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FIGURE  15.  Predicted  regions  of  positional  stability 
in  lift,  x/b=2. 
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FIGURE  17.  Predicted  regions  of  positional  stability 
in  pitch,  x/b=2. 
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FIGURE  16.  Measured  regions  of  positional  stability 
in  lift,  x/b=2. 


FIGURE  18.  Measured  regions  of  positional  stability 
in  pitch,  x/b=2. 


HASC95  predicts  four  distinct  regions  of  lift 
stability,  stable  in  lift  when  behind  and  below  the 
lead  aircraft  or  to  the  side  and  above.  These 
regions  are  divided  by  a  4-corner  boundary,  located 
at  0.40  span  spacing.  The  data  agree  well  with  the 
predicted  boundaries  with  two  small  exceptions. 
Instead  of  a  4-corner  boundary  between  regions,  the 
data  show  the  region  of  positive  stability  to  be 
continuous  with  a  small  channel  present  just  outside 
the  wing  tip  (0.6  span  spacing).  The  data  also  show 


Stability  boundaries  for  the  change  in  rolling 
moment  with  lateral  position  are  shown  in  Figures 
19  and  20.  HASC95  predicts  three  distinct  stability 
regions.  The  boundaries  of  these  regions  in  the 
horizontal  plane  of  symmetry  are  located  at  0.36 
and  0.94  span  spacing.  The  data  agree  very  well 
with  the  prediction. 
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FIGURE  19.  Predicted  regions  of  positional  stability 
in  roll,  x/b=2. 
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FIGURE  20.  Measured  regions  of  positional  stability 
in  roll,  x/b=2. 


In  all  three  cases,  there  is  a  break  in  the  stable 
regions  at  about  0.40  span  spacing.  For  the  basic 
horseshoe  vortex  model  of  a  wing,  this  is  the 
spacing  where  vortex  systems  of  the  wing  no  longer 
overlap  in  the  spanwise  direction,  i.e.,  the  trailing 
vortex  of  the  downstream  wing  moves  outboard  of 
the  trailing  vortex  of  the  wing  upstream. 

CONCLUSION 


Results  from  a  wind  tunnel  test  of  two  vertical 
tailless  delta  wing  aircraft  in  close  proximity  have 
been  presented  and  compared  with  predictions  from 
a  planar  vortex  lattice  method.  The  relative 
positions  of  the  aircraft  were  varied  in  the  spanwise 
and  vertical  directions  at  a  longitudinal  separation 
of  two  wing  spans  from  nose  to  nose.  Wake 
induced  effects  on  lift  were  large,  and  well 
predicted  except  when  the  aircraft  overlapped  in  the 
spanwise  direction.  In  that  case,  the  lift  increments 
were  overpredicted.  Wake  induced  effects  on 
pitching  and  rolling  moment  were  large,  and  well 
predicted  in  all  cases.  Wake  induced  effects  on 
side  force  and  yawing  moment  were  small, 
presumably  due  to  the  lack  of  a  vertical  tail.  The 
induced  drag  on  the  trail  aircraft  was  reduced  when 
the  aircraft  wing  tips  were  aligned  or  slightly 
overlapped.  Drag  increases  were  found  when  the 


wing  tips  overlapped  by  more  than  50%  span.  A 
maximum  induced  drag  reduction  of  25%  was 
measured  on  the  trail  aircraft  with  a  wing  tip 
overlap  of  15-20%  span.  This  compares  with  a 
predicted  reduction  of  40%  with  10%  overlap.  The 
discrepancy  was  attributed  to  flow  separation  at  the 
tip  caused  by  upwash. 
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APPENDIX  -  LONGITUDINAL  DYNAMIC  ANALYSIS 


Lor  formation  flight,  the  longitudinal  equations  of  motion  can  be  expanded  to  include  the  effects  of  vertical 
position  as  follows: 
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The  characteristic  equation  of  this  system  is  of  fifth  order  with  the  following  coefficients: 
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The  formation  terms  add  one  aperiodic  mode  to  the 
system.  Lor  stability,  a  necessary  (but  not 
sufficient)  condition  is  that  each  of  the  sn 
coefficients  be  positive.  A  quick  inspection  shows 
that  for  the  formation  derivatives  to  be  stabilizing, 
Lz  and  Me  should  be  negative  while  L0  should  be 
positive.  Mz,  which  only  appears  in  the  s°  and  s1 
coefficients,  presents  an  interesting  case. 
Intuitively,  one  might  expect  a  negative  value  is 
desired,  i.e.  the  aircraft  will  tend  to  pitch  towards  its 
original  flight  path  once  disturbed.  This  is  the  case 
for  the  s1  coefficient.  Lor  the  s°  term,  however,  a 
negative  value  is  destabilizing  for  the  following 
condition: 

cLcDa  >  cD(cLa  +CL9 ) 

Assuming  a  parabolic  drag  polar  and  CYe  =  0,  this 
condition  is  satisfied  if  the  lift  coefficient  is  larger 


than  that  corresponding  to  the  maximum  lift  to  drag 
ratio,  i.e.: 

CL  >  ARC Do 

Numerical  results  were  generated  for  a  full  scale 
version  of  the  tailless  aircraft  (b=37.5  ft)  using 
wind  tunnel  derived  stability  derivatives.  Pitch 
damping  was  computed  using  HASC95.  Two  flight 
conditions  were  examined,  M=0.6  at  10,000  and 
30,000  ft.  These  bracket  the  effect  of  Mz  on  the  s° 
term.  The  analysis  was  conducted  with  the 
formation  specific  derivatives  varied  one  at  a  time, 
at  their  maximum  values,  shown  below: 

CLz  =  ±0.0053, 1 1  ft  @y/b  =  0 
Cmz  =±0.0017,  lift  @y/b  =  0 
CLq  =  -0.138, 1  trad  @  ylb  =  0.8125 
Cmg  =0.046,  l/rad  @y/b  =  0.8125 
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Eigenvalues  for  sample  calculation,  M=0.6,  30,000  ft 


Clz 

cmz 

Clg 

Cm0 

A,  1,2  (short  period) 

A3, 4  (phugiod) 

A5 

-0.005 

-0.729  y  1.172i 

-0.012  y0.671i 

-0.013 

0 

0 

0 

0 

-0.744  yl.077i 

-0.004  y0.065i 

0 

0.005 

-0.711  yl.OOOi 

-0.779,  -0.013 

0.719 

-0.002 

-1.914  y  1.7181 

1.163  yl.636i 

0.005 

0 

0 

0 

0 

-0.744  yl.077i 

-0.004  y0.065i 

0 

0.002 
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0 

0 

-0.138 

0 

-0.732  yl.072i 

-0.015  y0.066i 

0 

0 

0 

0 

0.046 

-0.775,-1.690 

-0.002,  0.971 

0 

Eigenvalues  for  sample  calculation,  M=0.6,  10,000  ft 


CLz 

cmz 

Clg 

Cm0 

A  1,2  (short  period) 

A3.4  (phugiod) 

A5 

-0.005 

-1.555  yl.647i 
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-0.014 

0 

0 

0 

0 
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-0.007  y0.053i 

0 

0.005 

-1.606  y  1.517i 

-0.014,  -0.613 

0.655 

-0.002 

-2.712  y  2.198i 

1.127  y  1.948i 

-0.014 

0 

0 

0 

0 

-1.585  yl.584i 

-0.007  y0.053i 

0 

0.002 

-0.921  y  2.898i 

-3.331,-0.014 

2.002 

0 

0 

-0.138 

0 

-1.585  y  1.584i 

-0.007  y0.053i 

0 

0 

0 

0 

0.046 

-1.885  y0.958i 

-0.013,  0.599 

0 

The  results  are  shown  above.  The  height  stability  by  the  analysis  of  the  coefficients  of  the 
derivative  CLZ  has  no  effect  on  the  short  period  characteristic  equation.  None  of  the  Cm/  values 

characteristics.  It  does  affect  the  phugoid,  with  an  studied  result  in  a  completely  stable  system, 

unstable  value  changing  the  response  to  two 

aperiodic  modes,  both  of  which  are  stable.  The  The  effect  of  the  moment  derivative  Cm9  is  flight 

primary  effect  of  Clz  is  to  add  a  new  mode.  For  a  condition  dependent.  As  with  CLe,  it  almost  always 

positive  value  of  CLZ,  the  new  mode  is  unstable  appears  summed  with  the  pitch  stiffness,  Cma,  in  the 

with  a  time  to  double  amplitude  of  about  1  second.  characteristic  equation.  However,  it  is  twice  the 

magnitude  Cma,  so  it  has  a  significant  impact.  It 
The  lift  derivative  CL@  has  almost  no  effect  on  the  drives  the  phugiod  mode  to  two  aperiodic  roots,  one 

dynamic  modes.  It  almost  always  appears  summed  unstable  with  a  time  to  double  amplitude  of  about  1 

with  the  lift  curve  slope,  CL„,  in  the  characteristic  second.  For  the  low  altitude  flight  condition,  it 

equation.  It  is  only  6%  of  the  magnitude  of  the  lift  increases  the  short  period  damping,  and  for  the  high 

curve  slope,  however,  so  its  contribution  would  be  altitude  condition,  it  drives  the  short  period  to  two 

expected  to  be  very  small.  stable  aperiodic  roots. 

The  moment  derivative  C^z  has  a  major  effect  on  If  the  derivatives  were  varied  in  combination  with 

all  of  the  derivatives.  Although  the  short  period  each  other  and  not  one  at  a  time,  their  effects  on  the 

remains  stable  for  all  values  of  C^,  the  damping  is  dynamic  modes  would  probably  change, 
significantly  reduced  for  positive  C^.  Negative 
values  of  Cmz  cause  the  phugiod  to  go  unstable, 
with  a  time  to  double  amplitude  of  about  0.6  sec. 

Positive  values  drive  the  phugiod  root  to  an 
aperiodic  pair.  The  new  mode  is  unstable  for 
positive  Cmz,  with  a  time  to  double  amplitude  of  1/3 
sec.  The  effect  of  negative  Cmz  on  the  new  mode  is 
a  function  of  the  flight  condition,  as  hypothesized 
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